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1 Introduction 



The existence of plane-fronted parallel wave solutions in general relativity has 
been known since a long time ]j]]-pj]- These solutions are not only interesting 
by themselves, but also because similar solutions accompany massless fields of 
lower spin, such as electro-magnetic, scalar and spinor fields M, 
Gravitational pp-w&ves can be described by space-time metrics 

g^ v dx^dx v = — dudv — K(u, x, y)du 2 + dx 2 + dy 2 , (1) 

where the light-cone co-ordinates (u,v), transverse to the x-y-plane, are related 
to time- and longitudinal co-ordinates (t, z) by: 

u = ct — z, v = ct + z. (2) 

There are similar solutions with the roles of v and u interchanged, which we 
will not discuss explicitly here. With the metric (|l]), the only non- vanishing 
connection co-efficients are: 

p v — is y x_ in »_ \ K r S_ir "-lif (•}) 

1 uu ~ xv ,ut J- uu ~ 2 xu ~ 2 u ~ 2 yu ~ 2 

The free gravitational pp- waves are solutions of the free Einstein equations, 
which for the metrics ([!]) simplify to 

AtransK = K jXX + K >yy = 0. (4) 

The simplest non-trivial one is 

K(u,x,y) = ^k + (u)(x 2 -y 2 ) + K x (u)xy. (5) 

Here k^ + ^ x ^(u) are the amplitudes of the two different polarization states, as 
appropriate to quadrupole-type waves propagating at the speed of light. 

Metrics of the type (|l|) are also found in the solution of the coupled Maxwell- 
Einstein, massless Klein-Gordon-Einstein and massless Dirac-Einstein equa- 
tions. As a physically interesting example, consider the gravitational field of a 
light wave propagating along the z-axis. The light-wave is characterized com- 
pletely by a vector potential 1-form 

A = sin ku {a x dx + a y dy) , (6) 

with (a x ,a y ) the components of the constant transverse polarization vector. 
The corresponding electric and magnetic fields are 

E(u) = wa cos ku, B(-u) = k x a cos ku, (7) 

where oj = kc. These fields provide a solution of the coupled Maxwell-Einstein 
with a metric of the form ([!]), with 

K(u,x,y) = ^ (E 2 (u)+c 2 B 2 (u)) (x 2 + y< 



1 



The motion of a particle of mass m and charge q propagating in gravitational 
and electro-magnetic background fields is governed, in the limit of neglecting 
radiation reaction terms, by the covariant Lorentz force 

x» + T x /± x x u = F%x v . (9) 

Here the overdot denotes a derivative w.r.t. proper time, defined in the usual 
way from the universal conserved value of the total four-velocity: 



2 dap dxt 1 2 

U » = 9 ^^=- C - (10) 



We apply these equations to study the motion of a particle in the background of 
the waves fields (||J6|,0) • Recalling the results of ref. [0 , the light-cone co-ordinate 
u can be used as a proper-time variable via the relation 

u = 7 = const. =>■ u = 77". (11) 

The expression ([!(]) for the total four-velocity then leads to the following rela- 
tion between proper time and laboratory time: 



h = K +T, Z J T^ = t 2- ( 13 ) 



dt _ dt __ l-~i 2 K/c 2 

~ ^ ~ V 1 " v 2 /c 2 ' ( } 

showing that there is both a gravitational and a kinematic time dilation. Sub- 
stitution back into eq.([To|) allows one to cast it into a true conservation law of 
the form 

1 _ v 2/ c 2 _ c 2 

(1 - v z /cf ~ 7 2 

Two independent equations of motion remain, which we take to be those for the 
tranverse co-ordinates £ = (x,y). The Lorentz-Einstein equation (0) for these 
components becomes 

—4 = 7T~ k 2 a? (1 + cos2£;u)£ - — kacosku. (14) 

cm 2 cr 7777 

We can define the components w.r.t. the polarization of the light wave as 

£ = £u + £±, Z\\ = rw a > Z± = ^rr- ( 15 ) 

1 11 |a| z |a| 

Under this decomposition eq.(|l4]) splits into a homogeneous and an inhomoge- 
neous equation 

+ (1 + cos 2ku) € ± = 0, 

" J G (16) 

J -1 ^— /c 2 a 2 (1 + cos 2ku) £ii = — — fcacos ku. 

du z c z " 7717 
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In refs. |], [?J it was observed, that in the static limit k — > taken such that 
Eq = cBq = ck\a\ = constant, one obtains a simple harmonic motion with 
frequency v (in Hz) given by: 



2-kv = W J E = 0.3 x KT 18 £ (V/m). (17) 
In this paper we discuss the motion in the non-static case k ^ 0. 



2 Transverse motion 



We first focus on the case q = 0, describing a neutral massive test particle in an 
electro-magnetic and gravitational pp-wave. According to (|l6|), the transverse 
coordinates (x, y) all obey the same homogeneous Mathieu equation. Besides, 
we choose to restrict our study to the case where the only contribution to the 
metric co-efficient K is the electromagnetic one, given by Iffi). In other words, 
no free gravitational pp-waves are taken to be present: k+ * (u)=0. It then 
follows that the motion of the particle is invariant under the exchange of the 
tranverse coordinates and rotations in the transverse plane. The equation for 
£ = {x,y) now reads (with a 2 = a 2 ): 



d 2 i 27T£og ,22/-,, , w n 

-t-o H o — k a (1 + cos 2/cu) t = 0, 



(18) 



A solution procedure for eq.(18) in the weak-coupling regime was presented 
in ref. ||. Here we discuss the construction of the full solution in terms of 
a generalized Fourier series expansion. Because of the particular form of the 
equation (^) , we can construct a series in terms of either even or odd multiples 
of the fundamental wave factor k. In this section we consider the generalized 
odd solutions, which can be parametrized as 



oo 

= e aku (£n cos(2n + l)ku + rj n sin(2n + l)ku) , 

n=0 



(19) 



with real Fourier components (£ n , rj n ). Note that for cr 2 < (and taking the real 
part) one obtains solutions which are purely of trigonometric type; for a 2 > 
the solutions describe oscillations enhanced by parametric resonance. The other 
type of solutions having even Fourier part is discussed in the appendix. The 
substitution ([l90 converts the Mathieu equation into an infinite-dimensional 
linear matrix equation 



/ A 
B 


V • 



B 
B 




B 
A 2 





( Po\ 




pi 




P2 


) 


{ ■ J 



(20) 
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where the blocks A m and B are 2x2 matrices given by 



A m 

and 



a- 



1 + 
-2a 



3k^ 



<J- 



2a 

1 ~ 4 



(21} 



a 2 - (2m + 1) 2 + f 



2(2m + l)a 



-2(2m + l)<r a 2 - (2m + l) 2 + \ 



(m > 1), 



re 

T 



1 o 
1 



(22) 



A solution of the system ( J20| ) exists provided the linear matrix operator has zero- 
modes. To determine their existence we cast the equation into upper-triangular 
form: 



ax B 
a 2 



V 





( Po 


\ 




Pi 






P2 




) 


K ■ 


J 



0. 



(23) 



The 2x2 matrices be expressed in terms of continuing fractions as 

-l 



Am - B z 



B 2 [A m+2 



and the basis of vectors has been transformed by po = POi an d 
pm = Pm + a^Bp m -i, (m > 1). 



(24) 



(25) 



This transformation is singular if and only if any of the a m (m > 1) has zero- 
modes, which is precisely the condition for the Mathieu equation to have non- 
trivial solutions. 

Zero-modes of ao represent the lowest-frequency solutions of the Mathieu 
equation; they exhibit parametric resonance for a > 0. The existence of zer- 
modes implies that the determinant of ao vanishes. This determinant can be 
computed in successive approximations by truncating the infinite fraction (E4|) . 
The zeroth-order approximation is 



ao = A 



,(°) 



{a 2 + l) 2 + K 2 (a 2 -l) + 



3k 4 
Iff 



The first- and higher-order approximations are of the form 



a 



l + 3|! + /t 4 An 



2_1_ L 3«"_ 1 ..4 a.. 2<T + K 4 Ai2 



■(2a + K 4 A 12 ) a 2 



1 + \ + k An 



+ C(k 6 



(26) 



(27) 



where the quantities Ay(a) are rational functions of a 2 ; the determinant then 
is 



Nl 



(a 2 + l) 2 + K 2 (a 2 -l) + K 4 



A 



16 



+ 0(k 6 ). 



(28) 
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In first order = A$ — B 2 A 1 1 , with 

(1) _ l_ 9-a 2 (1) _ 1 1 

11 48 27 + 6a 2 + a 4 /3' ^ 12 8 27 + 6a 2 + a 4 /3 ' 1 j 

and 

A (D - ^ 4 - 22 ^ + 9 (30) 
Similarly, in second order with 

a = (AqA^ - B 2 (A + A 2 j) (A^ 2 - B 2 ) ~* , (31) 

one finds up to terms 0(k 6 ): 

/9 , A ( 2 ) 5 3 (a 2 + 25)(a 4 -94a 2 + 225) 



(cr 4 - 94a 2 + 225) 2 + 256a 2 (a 2 - 15) 5 

A (2) 3 (a 2 + 25) 2 

12 8 (a 4 - 94a 2 + 225) 2 + 256a 2 (a 2 - 15) 2 ' 



(32) 



Then in the determinant 

A(2 ) = (a 2 - l)[(a 2 - 25) (a 4 - 94a 2 + 225) + 160a 2 (a 2 - 15)] 

(a 4 - 94a 2 + 225) 2 + 256a 2 (a 2 - 15) 2 ' 1 j 

Observe that in both first and second order approximation 

lim A^ = -, lim A^ = 1, k = (1,2). (34) 

a 2 — >0 9 <r 2 -^oo 

From the above analysis it follows, that up to terms of order k 6 a solution for a 2 
exists for which the determinant vanishes: |ao| =0, as given by the expression 
(f2Sl). This condition has approximate solutions 



K 2 



a 2 = -[!+■-) ± ^ 2 (l + ^(l + 2A)). (35) 



2/t 2 I 1 + !L (1 + 2A) I > I 1 + '— I . (36) 



Positive non-zero values of a exist only for the positive square root, provided 

3~2 
or 

fl-y +^VYT2a)(i-^-^VTT2a)<0. (37) 



Hence k 2 must be in the range 



1 K, 2 1 

< — < . (38) 



2 + V1 + 2A 4 2- V1 + 2A 
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Purely periodic solutions of the Mathieu equation exist for a discrete set of 
values of k for which a = 0; for these values the inequality (^) turns into and 



equality. For all other values outside the range (38) the solutions for a are 



imaginary; the series (19) then is purely trigonometric, but for non-rational a 
the result is not periodic. 

From the inequality ( |38[ ) it follows, that the smallest possible value of k 2 
for which parametric resonance occurs is of the order unity: k 2 > 1. As k 2 
represents the ratio of field intensity and frequency, it follows that one needs 
very strong fields at very low frequencies (2ttv = ck): 

f0.3 X IQ-^o (V/m)\' 

^ 2vr^(Hz) J ~ v ' 

An upper limit to the fields strength of the order of 10 19 V/m arises because 
of pair creation processes which destroy the electric field at higher intenties |§ . 
Therefore we obtain an upper limit of the order of ~ 1 Hz for the frequencies 
which can give rise to parametric resonance. For lower fields and higher fre- 
quencies one finds standard periodic behavior described by solutions with a = 
or imaginary. 



An important characteristic of the motion of a chargeless particle in the back- 
ground of the fields (|l],^,@) is that a particle initially at the origin in the trans- 
verse plane before the arrival of the wave, remains at the origin at all later 
times. This allows one to translate directly the geodesic motion of a test par- 
ticle in terms of the relative rate of acceleration of this particle w.r.t. another 
one located in the origin: £ = 0. 



For charged particles this situation changes. If q ^ and we take the transverse 
co-ordinate £ = £y parallel to the electric field, the equation of transverse motion 

1 ■•>'■>' -e cos ku, 



dv? 



+ 



(1 + cos 2 Ant) £ 



with 



qka 

rwy 



mcj 



(40) 



(41) 



In the inhomogeneous case the absence of an expontential term in the driving 
force implies that any special solution of the inhomogeneous equation can have 
an expansion ( |i~9D only with a = 0; for these special solutions the inhomoge- 
neous Mathieu equation becomes 






V • 



B 

at 





B 

(12 



{ Po\ 

pi 

P2 
\ ■ J 



( o \ 




V • / 



£ 



l 

o 



(42) 



Here we are interested especially in obtaining a particular solution of the pre- 
vious type: po ^ 0, p n = 0, n > 1 . The non-zero component po is a solution of 
the inhomogeneous equation 

a p = 0. (43) 
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Now ao is of the form (pTj), with a = 0: 



«o 



l + ^ + A n (0)K 4 





+ £ + A n (0) K 4 



(44) 



where beyond the lowest-order approximation we find An (a = 0) = tjt. As a 
result a special solution of the inhomogeneous equation is given by 



P0 



k 2 (l 



3i£ 
4 



Aii(0)k 4 + ...) 



1 





0, n > 1. 



(45) 



The general solution of this system is the sum of a particular solution plus a 
solution of the homogeneous equation, which can have a / 0. An important 
difference with the previous case is that a particle initially at the origin does 
not remain at the origin, but is accelerated by the electric field. Thus the study 
of the relative acceleration between particles on different world lines becomes 
more complicated; for a general discussion we refer to ref. Jul. 



3 QED corrections 

As the magnitude of the electric field in the region of parametric resonance 
is necessarily extremely high, QED effects such as vacuum polarization can in 
principle not be ignored. However, as we now show they do not affect the 
classical wave solutions we have studied in the previous paragraphs. 

The non-linear effects of quantum fluctuations on the propagation of light 
in flat (Minkowski) space-time were first studied by Euler and Heisenberg ||. 
The effects are summerized by a well-known effective lagrangean incorporating 
the corrections due to the electromagnetic higher-order terms; its generally 
covariant form reads: 



:„c\/—[~F 2 + (3 



where: 



g^g vX F. v F KX , 



(46) 



(47) 



and Fn U is the dual electromagnetic tensor defined by: 



F» v = e^^Fpv. (48) 
The coupling constant f3 is given by 

P = ^f^. (49) 

H 90m 4 c 4 v ; 

We first derive the Maxwell equations from the Euler-Lagrange equations taking 
A v and d u A^ as the variables: 

if-^TO^ - (5o) 
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Combining (46) and (pOl), we get the following equations: 



V, 



(1 - 8/3F 2 )F^ - Uf3<pF 



0. 



The energy-momentum tensor is obtained from the lagrangean 
-2 



1 fJ,U 



-9 dg 



e c 



by: 



(1 - 8(3F 2 )F, X F U X - i^(F 2 - 4/?(F 2 ) 2 ) - ^ g^ 



(51) 



(52) 



With the metric @) the Ricci tensor R^ u has only one non-zero component R uu , 
the Einstein field equations reduce to: 



1 



AtransK = 8irGT u 



all other = 0. 



(53) 



In our particular case of a plane- wave vector potential of the form (^), the only 
non-vanishing components of the electromagnetic tensor F„ v are: 



F„ 



-F 

± 77 



dik cos ku, 



(x,y)- 



(54) 



It follows that F 2 
the original form: 



0, and the Maxwell-Einstein field equations reduce to 



V M F^ = 0, 



2 A 

-Eqc F u \F u . 



(55) 



Thus the quantum corrections do not modify the gravito-electro-magnetic plane 
wave solution. 
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Appendix 



The solutions of the Mathieu equation are classified as having an odd or even 
periodic part. Eq.(|l9|) represents the general solution of odd type. The general 
solution of even type can be expanded as 



£(u) 



Mo 
2 



oo \ 

+ ^ [fi n cos 2nku + v n sin 2nku] ) . 



n=l 



(56) 



After multiplication with e uku the homogeneous Mathieu equation becomes 

,.2 







1 2 1 2 

- cr +-k I Mo + ~r Mi 



a 2 - 4n 2 + — I /x n + 4no-^„ + '— (fj, n -i + Mn+l] 



+ E 

n=l 

oo 

+ E 

n=l 

+ ( ^2 sin 2ku + ^ (z^i-l + v n+l) si n 2n/cu 



cos 2n/cu 



a — An H — — I v n — Anafi n 



sin 2nku 



It follows, that 



and 



n=2 



Mo 



^Ml 

2a 2 + k 2 ' 



(D 1 


5 









\ 


B 


D 2 


B 




Ci 







D 


D 3 




C2 




V • 






• / 


V • 





(57) 
(58) 

(59) 



where B is the 2 x 2-matrix (p2[), whilst 



L>1 



A, 



a 2 - 4 + £ 



4(2o- 2 +k 2 ) 



and 



-4<T 



o- 2 - 4n 2 + ^ 
— 4n<r 



4cj 

a 2 - 4 + £ 



(60) 



4n<7 
a 2 - An 2 + ^ 



2 



n > 2, 



Mn 



(61) 



To solve this equation, we bring the infinite-dimensional matrix equation (|59|) 
into upper triangular form: 



( D x B ... \ 
B D 2 B ... 
B D 3 ... 



\ 



J 



( d x 


B 








d 2 


B 








d 3 


V • 







( 1 ... \ 
ci 1 ... 
c 2 1 ... 



V 



(62) 
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with d n and c n the solutions of the equations 

dn+Bcn = D n , dn+icn = B => B^d n + {B~ l d n+l )- 1 = B~ l D n . (63) 
The solution for d n takes the form of a continuing fraction 



d n = D n — B 



D n+1 - B 2 \D n 



+2 



(64) 



where we have used the fact that B is proportional to the unit matrix, hence 
commutes with all other matrices. If we define 



(2 
(3 



( 1 

a l o 

c 2 1 



Ca 



(65) 



• / 


V • 




• / 


\ ■ J 


, and 










Cn = 




lCn-l 


n 


> 2. 


/ di 


B 


... \ 


(^\ 







d 2 B 




C2 


= 0, 





d 3 






V • 




• / 


V • J 





Now the equation 



(66) 



(67) 



has solutions with (jv 7^ 0, and 

d k (k = -B(k+l, l<k<N-l and (n = 0, n > N. (68) 
The first one of these is the one given by di(i = 0, ( n = 0, n > 2. 

Approximations. Solutions of the equation di^i = exist if det d\ = 0. To solve 
this condition, we need to compute the infinite fraction (|64|). We can do this 
by successive approximations: 



d\ 



(0) 



7(1) 



D x - B 2 D 2 l = (DiD 2 - B 2 ) D 2 \ 



4 2) = (D^Dz - B 2 D 1 - B 2 D 3 ) (D 2 D 3 - B< 



(69) 



To zeroth order approximation, we can show that det D\ = has solutions for 
8 < k 2 < 16: 

\Di\ = y( k2 - 8 ) (* 2 " 16 ) + 2a2 ( a4 + ( 8 + ¥) ° 2 + 16 + • (?0) 

It suffices to note, that the first term is negative in the range k 2 £ (8, 16), whilst 
for any k 2 the second term can take all positive values between (0, oo). 



10 



References 

[1] M.W. Brinkmann, Proc. Natl. Acad. Sci. 9 (1923), 1; 
Math. Ann. 94 (1925), 119 

[2] H. Bondi, Nature 179 (1957), 1072 

[3] J. Ehlers and W. Kundt, in: Gravitation: an Introduction to Current Re- 
search (Wiley, N.Y.; 1962), ed. L. Witten 

[4] D. Kramer, H. Stephani, M. MacCallum and E. Herlt, Exact Solutions of 
Einstein's Field Equations (Cambridge Univ. Press; 1980) 

[5] J.W. van Holten, Phys. Lett. B352 (1995), 220 

[6] J.W. van Holten, Fortschr. Phys. 45 (1997), 6 

[7] J.W. van Holten, in: Towards Quantum Gravity, Proc. XXXV Int. Win- 
terschool on Theoretical Physics (Polonica, 1999); (Springer, 2000), 365 

[8] J. Schwinger, Phys. Rev. 82 (1951), 664 

[9] W. Heisenberg and H. Euler, Z. Physik 98 (1936), 714 

[10] A. Balakin, J.W. van Holten and R. Kerner, Class. Quantum Grav. 17 
(2000), 5009 



11 



